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Abstract 

Overlapping resonances are shown to provide new insights into the extent of decoherence expe- 
rienced by a system superposition state in the regime of strong system- environment coupling. As 
an example of this general approach, a generic system comprising spin-half particles interacting 
with a thermalized oscillator environment is considered. We find that (a) amongst the collection 
of parametrized Hamiltonians, the larger the overlapping resonances contribution, the greater the 
maximum possible purity, and (b) for a fixed Hamiltonian, the larger the overlapping resonances 
contribution, the larger the range of possible values of the purity as one varies the phases in 
the system superposition states. Systems displaying decoherence free subspaces show that largest 
overlapping resonances contribution. 

PACS numbers: 03.65.Yz,03.67.Pp 
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I. INTRODUCTION 



Coherent control of atomic and molecular processes is a subject of considerable ongoing 
interest . Specifically, a number of successful control scenarios have been theoretically 
proposed and experimentally implemented. The vast majority of well characterized scenarios 
are, however, based upon molecules treated in isolation, with the understanding that systems 
in an external environment are subject to decoherence effects[3,kfl that can be deleterious to 

n 

control|5[. Decoherence effects of this kind are also detrimental to "quantum technologies", 
such as quantum information processing and quantum computingjg] insofar as they cause 
the loss of the essential quantum characteristics that make these research areas interesting. 

As a consequence, efforts are ongoing to find suitable mechanisms to reduce or compensate 
for decoherence effects. For example, several techniques, such as quantum dynamical de- 
coupling (QDD) [7j and the use of decoherence free subspace (DFS) [8], have been proposed 
in the quantum information community. However, the QDD method is quite challenging to 
implement using the present technology, and the DFS requires explicit symmetries in the 
Hamiltonian, which are not generic. Hence, methods for identifying conditions for states 
that are resistant to decoherence is of great interest. 

Indeed, considerable effort has gone into identifying characteristic features of particular 
system states, interacting with different environments, that lead to stability against deco- 

lerence. Sample approaches include the general predictability sieve approach of Zurek et al. 

9] , and predictions of the unusual stability of coherent matter states, low lying eigenstates, 



and localized states 
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11] for typical environmental coupling. 



In this paper we introduce a new perspective, based on overlapping resonances, to under- 
standing and optimize the resistance of system superposition states to decoherence. The ap- 
proach applies in the important regime where system-environment coupling is large. Specifi- 
cally, we show that interference between the overlapping resonances that exist in this regime 
can be used to build system superposition states that are relatively resistant to decoherence. 
The extent to which such states are stable against decoherence is shown to directly correlate 
with the degree to which the system-environment Hamiltonian displays overlapping reso- 
nances. Thus, even if the system does not possess symmetry properties such as decoherence 
free subspaces (DFS), states can be constructed which are more resistant to decoherence 
than are others. 
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From a formal viewpoint, the approach described herein constitutes a major extension, 
to product Hilbert spaces characteristic of open quantum systems, of our work on control- 



ling radiationless transitions in large molecules (12J. There, population transfer between 
a Q-subspace and the remainder of the Hilbert space, the P-subspace, was minimized by 
preparing initial superpositions in the Q-subspace, relying on the presence of overlapping 
resonances. This approach was applicable to the case where the full Hilbert space was par- 
titioned into the sum of two subspaces. By contrast, decoherence involves the effect of an 
environment on the system, i.e. a product of two Hilbert spaces. The difference is substantial 
and significant. 

The inter-relationship between overlapping resonances and the ability to design states 
that are resistant to decoherence is general. It proves useful, however, to introduce the 
approach via a specific generic example. For this reason we choose a particular class of spin- 
boson problems. Note that system decoherence is a manifestation of the loss of quantum 
information when one neglects the environment with which the system is entangled)^]. As 



such, decoherence can occur effectively even with a small environment 13J] as is manifest in 
the example below. 

The structure of the paper is as follows. In Sec. II, we introduce a generic system that 
describes a spin system interacting with a environment, providing a specific example for 
this general approach. Section III, considers the nature of decoherence in this system, and 
its relationship to overlapping resonance. Some specific numerical examples are provided 
in Sect. IV to demonstrate correlations between overlapping resonances and reduced state 
decoherence. Sect. V provides a summary. 

II. A GENERIC CLASS OF HAMILTONIANS 

The spin-boson problem has been extensively studied in Chemistry in various scenarios 



and for a wide variety of applications 14j. Here we consider, as a specific example, the control 
of decoherence in a simple spin system interacting with a small thermal environment, i.e. a 
single spin-half system and an environment comprising a single-mode oscillator. The most 
general system-environment Hamiltonian for this case is given by 

H/h = ujo)a + ~^S Z + g r (S + a + S^a^) + g nr (S + a) + SLa) + g p hS z (a + a)) , (1) 
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where a and a' are the environment annihilation and creation operators, u is the frequency 
of the environment mode, hwo is the energy difference between the two spin-half states |±), 
and the spin operators are defined as S± = |±)(=p|, S z = |+)(+| — |— )(— | • Note that the 
rotating wave approximation (RWA) has not been invoked: the g r term in the Hamiltonian 
is the spin-environment interaction that is retained under the RWA, whereas the g nr term is 
that which would be neglected in the RWA. The g p h term corresponds to decoherence of the 
spin due to the operator S z , which induces a relative phase between the spin states |+) and 
|— ). The g r and g nr terms mediate energy exchange between the spin and the environment. 
Our focus below is in the region of strong coupling, where g nr and g r are both greater than 
ojq. It is important to emphasize that, although the environment is small (a single oscillator) 
its coupling to the system does induce significant system decoherence. 

The class of Hamiltonians [Eq. <^j] possesses a plane defined by g r = g nr in the three 
dimensional g n r,gr,g P h parameter space upon which the system displays a decoherence free 
subspace (DFS). That is, in these special cases the Hamiltonian Eq. (CQ) leads to a DFS due 
to a Lie algebraic symmetry. In these cases, the interaction part of the total Hamiltonian 
can be written as Ei ® Bj, where Bj = (a + a 1 ) is the environment operator, and the spin 
operator Sj can be written in a closed matrix form in the (|+), |— )) spin basis: 



\ 9 -9 P h J 

The DFS corresponds to the eigenstates of the Si matrix. For other values of g nr and g r , of 
primary interest here, the system does not have a DFS. 

III. DECOHERENCE OF A SPIN SUPERPOSITION 

We consider the stability against decoherence of an initial spin superposition state 
l^) = c + |+) + c_|— ), where c± are complex amplitudes and the environment is in thermal 
equilibrium at temperature T. The latter is described in terms of the oscillator eigenstates 
|n) by the environment density matrix ps = Yl n Pn\ n )( n \: Vn — e~ En l kBT e~ En ^ kBT ) and 
E n = (n + 1/2) tvjj. In order to follow the temporal evolution of the system we expand the 
total (system+environment) evolution operator U in terms of I7), the eigenstates of the full 
Hamiltonian [i.e., (E^ — H)\^) = 0] with U = e ~ lE i t / h ^ Expanding the eigenstates 
1 7) in the zeroth order basis of "bare" spin+environment states {|±, n) = |±) <g) (i.e. 




(2) 
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the eigenstates of H without the spin-environment interaction) gives, for the time-evolved 
density matrix of the total spin+environment system: 

p(t) = 5> n f/[|^| ® |n><n|] U^=J2 Pne^-^h)W\(^,n)(^,n\i) • (3) 

n 7,7',n 

Further, tracing over the environment states, gives the following reduced system density 
matrix elements (p s ) k i in the basis of spin states \k) and \l) where k, I denote either the plus 
or minus state: 

(Ps)m= Pne-^- £ v)*(fc,m| 7 )(7^,n)(^n| 7 ')(7'|Z,m) . (4) 

m,n,7,7' 

Expanding as c + |+) + c_|— }, we can rewrite the above expression as 

(p s ) kl = \c + \ 2 Q kl {t) + \c„\ 2 R kl {t) + c+c*_P kl (t) + c^c* + T kl {t) , (5) 
where P,Q,R,T are complex matrices with elements 



Pklit) 
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Consider the decoherence of the spin state, quantified in terms of the purity 
spin, where S = Tr(p 2 ). Given the above results, S can be evaluated via 



(6) 

S of the 



s= E (*wxi|p.i*> = E i(*)«i 2 ( 7 ) 

fc,ie+,— fc,ze+,- 



so that, using P^ = X^: 

s = | c+ i 4 ]T |g M | 2 + |c-i 4 ]T |p* 

k,le.+ - k,ie+,- 
+ |c + | 2 |c_| 2 £ [|P w | 2 + |T w | 2 + {Q fc ^ + Q*^}] 



+2Re 



£ [c 2 + cI 2 P M P fc + \c + \ 2 clc4Q kl P^ + P^Ql,} + M^c^P* + P; fc P* z ]] (8 



.fc,/6+, 



The first three of these terms depend solely on the magnitude of the c± coefficients, whereas 
the remainder of the terms depend upon the phases of the coefficients as well. The latter 
are clearly contributions to S that rely upon quantum interference between the spin states. 

Note, for comparison later below, that S — 1 denotes a pure state. If the state is 
completely decohered in the spin basis then (p s )ki = \&ki so that the purity, by Eq. ([7]), is 
S = 2(l/2) 2 = 1/2. 

With the P, Q, R, T matrices known, one can optimize S with respect to c + and c_ to 
determine the state with maximum purity S at a fixed time. In particular, the relative phase 
of c + and c_ is important if either Qm or Rki is non-zero, and Pm is non-zero . Although not 
evident from Eq. ([S]), numerical results below clearly demonstrate a reliance on this relative 
phase to maximize S. 

IV. OVERLAPPING RESONANCES AND SYSTEM PURITY 
A. Overlapping Resonances 

As noted above, the bare states |±,n) are not eigenstates of H, and will evolve in the 
presence of the system-environment coupling. Insight into the nature of this time evolution 
is afforded by expanding the zeroth order state in the exact eigenstates I7). For example, a 
bare state such as |+, j) expands as 

I+,.7> = £|7><7|+,J>, W 

7 

and evolves as 

7 

The E 1 dependence of the square of the expansion coefficient, i.e. |-D 7 | 2 = |(7|+,j)| 2 
provides the energy width over which the zeroth order state \+,j) is spread due to the 
system-environment interaction. In the simplest cases the inverse of this width provides a 
qualitative measure of the time scale for the evolution of |+, j). Hence, the zero order states 
are indeed resonances with a characteristic width |-D 7 | 2 . This is the analog, in a bound state 
spectrum, of the well known resonance in the continuum. 

Similarly, by analogy to the continuum case, we can define overlapping resonances, as 
resonances that overlap in energy space, i.e. resonances that share a common I7) in their 
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respective eigenstate expansions. These are the precise bound state analogues to overlapping 
resonances in the continuum [e.g. Fig 1, Ref. 16|]]. Examples of bound state overlapping 
resonances in molecular systems have also been previously reported 12]. 

In the case of the spin-boson Hamiltonians given by Eq. ([I]), overlapping resonances 
were clearly evident in, but not limited to, the strong coupling regime characterized by g r 
and g nr greater than uq. They occur as two states \k,m) and \l,n) with k ^ I and/or 
m ytz n (k, I G +, — ) that have non-zero overlap with the same eigenstate I7). Examples 
of such overlapping resonances are shown in Fig. [TJ Quite clearly the overlap between the 
resonances can be significant, enhanced by the fact that the overall width of each resonances 
can be significantly larger than the energy spacing between the zeroth order states. The 
resonances are also seen to be rather highly structured. 



B. Contributions to Tr p 2 s 

Consider then the role of these overlapping resonances in the dynamical evolution of p s (t) 
given in Eq. (j3J). Specifically, note that all the matrices P,Q,R,T depend on overlapping 
resonance contributions. Their dependences differ, however, in the important case where 
k 7^ l,k,l G +,— • In that case overlapping resonances contribute to and P^i, but not to 
Qki and R^i- That is, the terms in (p s )ki that are sensitive to the relative phase of c + and c_ 
[see Eq. (jSJ)], and are hence associated with phase control over the dynamics, are directly 
dependent on overlapping resonances contributions. 

^From Eq. (jSJ), it is apparent that the purity is strongly affected by overlapping reso- 
nances through the matrix elements Qki, Rki, Pkh an d T^. To quantify this dependence, we 
explicitly identify overlapping resonances contributions associated with the plus and minus 
spin states. That is, we define 



A+ _ = J2 A rn,n A rn,n ^ ^.n 



^|<±,m| 7 )( 7 |±,n)| 



(11) 



For a pair of the oscillator states \m) and \n), A™ ,n determines the overlap (summed over 
all eigenstates I7)) between the bare states \+,m) and \+,n), when the spin state is |+). 
If, for the same pair of |m) and |n), there exists non-zero overlap A™ ,n between the bare 
state \—,m) and |— ,n), then A + _ provides a measure of the overlap between states |+) and 
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|— ), that arises through the overlap between states |±,m) and \±,n). Note that Eq. (|TT|) 

is somewhat similar to the off diagonal T + _ and P H terms at t — 0. 

Note that the number and character of the overlapping resonances are functions of the 
parameters in the total Hamiltonian, e.g., here g r , g nr , g p h, 00, and Uq. For example, in the 
extreme case where there is no system-environment coupling (i.e. g r = g nr = g p h = 0), the 
total Hamiltonian is diagonal in the spin-boson bare state basis. The overlap between the 
spin-boson bare states then vanishes, and in turn, so does the overlap between the states 
|±). That is, as expected, in absence the system of coupling to the environment, the spin 
states are not broadened and thus do not overlap. 

C. Computational Results 

Consider then the relationship between overlapping resonances as reflected in v4 + _ and 
decoherence as embodied in the system purity S. As an example, we examine dynamics 
where the environment temperature T = 25 111K and choose to vary g r , keeping the other 
Hamiltonian parameters fixed. 

Numerical results are obtained by diagonalizing the Hamiltonian matrix in the zeroth 
order basis of bare states |±,n). At the temperature of T = 25 mK, a total of 20 oscillator 
states contribute, allowing us to diagonalize a 40 x 40 dimensional Hamiltonian matrix to 
obtain the desired E 7 , 17). 

To assess the decoherence we find, for each value of g r , the maximum possible purity 
'S'max in the system at a fixed time t by optimizing the complex coefficients c±. Results are 
shown in Fig. [2] where both S majX and the overlapping resonance measure A + _ are plotted 
as a function of g r at t — 0.1 ns. The correspondence between the behavior of these two 
functions is evident, with larger A + _ correlating well with larger S m&x . That is, the greater 
the overlapping resonances contribution, the more resistant is the optimal superposition to 
decoherence. Note, significantly, that the optimally resistant state need not be, and is often 
not, an eigenstate of the spin Hamiltonian, but is rather a superposition of |+) and |— ) 
eigenstates. 

Prominent in the graph (and some presented below) is the DFS point g r = g nr , where 
the A + _ curve displays a sharp peak, commensurate with S mSiX reaching its maximum value 
of unity. That is, maximal purity clearly correlates with the maximum The deviation 



S 



of S max from unity as one moves away from this point is significant, particularly in scenarios 
such as quantum computation, where extremely high degrees of coherence are necessary. 

The results clearly show that states take advantage of overlapping resonances (and hence 
of quantum interference) to increase the state purity This is the case even if the system- 
environment coupling is very large, as is clear from Fig. [2(b), where, for example, g r ~ 
27r x 1.5 GHz (g r 3> ujq), but where the purity remains as large as 0.92 for g p h = 2n GHz. 

The presence of overlapping resonances, as seen above, influences the degree of decoher- 
ence when preparing initial states with different initial coefficients c±. Sample dependences 
of 5* on c± are shown in Fig. [3] for the case of g r = 2tt x 400 MHz, g nr = 2tt GHz, and 
g p h = 2tt x 500 MHz [i.e. parameters associated with Fig. [2(a)]. Shown are the time 
evolution of 5* for three states, one being the best case superposition c + = —0.539998 and 
c_ = a/ 1 — | c + 1 2 which shows S max = 0.96 at the time t = 0.1 ns where 5* is optimized. 
Comparison is made with the time evolution of S for two other sample states with different 
c± values: c + = 1 and c + = c_ = 1/ y/2. The range of S as a function of initial state is large, 
with the c + = 1 case showing almost fully mixed behavior at long time. Interestingly, the 
curves shown do not cross as a function of time. Hence the optimal superposition, deter- 
mined by coefficients obtained for the t = 0.1 ns case is, in fact, the optimal superposition 
for the times preceding that target time as well. This was also confirmed by optimizations 
independently carried out at shorter times. 

Additional results showing how overlapping resonances between spin states is manifest in 
decoherence is shown in Fig. [H Here, for a given value of g r , we obtain the \c + \ and |c_| that 
maximize S and then vary the relative phase 9 between c + and c_ to obtain Sdis, defined 
as the difference between the maximum and minimum purity so attained. This difference is 
plotted in Fig. H]as a function of g r , along with the corresponding value of the overlapping 
resonance contribution A + _. The extent to which Sdis varies is seen to reflect the variation 

of A + _. Specifically, as A] increases, so does the extent to which the purity varies with 9. 

Hence, in the presence of large overlapping resonance, one can actively control the purity of 
the system by taking advantage of the phase dependent quantum interference contribution 
to S. Indeed, in accord with Eq. ([8]), an observed variation of the decoherence with the 
phase of a superposition state provides evidence of the presence of overlapping resonances. 

The extent to which S^is varies with A + _ depends upon the Hamiltonian parameters. 
Consider, for example, an interesting case of the nanoscale Hamiltonian system in Eq. (CQ), 
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the interaction of a Cooper pair box (CPB) with a nanomechanical resonator [l 
system is given by: 



This 



Hcpb = 4:E c 5nS z - (Ej/2)S X + fkoa^a + hgS z (a + a f ) , (12) 

where Eq and Ej are the charging energy and the Josephson energy of the CPB, u is 
the fundamental frequency of the resonator, and Sn lies between —1/2 and 1/2. At the 
degeneracy point 5n = this Hamiltonian can be written in the form of that of Eq. (JT]) 
under a similarity transformation, where we identify uq = Ej, g p h = 0, and g r = g nr = g. 
The similarity transformation is given by the operator e », with S y being the y-projection 
of the spin operator. 

^From the symmetry of the Hamiltonian (TT2|) . one finds a DFS corresponding to the state 
c + = c_ = 1/V2 and c + = — c_ = l/y/2. Hence, the value of S m3iX will always be unity for 
all values of g, irrespective of the values of A + _. Nonetheless, overlapping resonances affect 
the decoherence in this system. Specifically, as A + _ increases with g, the variation of S with 
changes in the relative phase 9 between c± also increases. This is shown in Fig. [5] where the 

difference S^is between the maximum and minimum values of S and the values of with 

g are seen to be strongly correlated. 

On the other hand, if g r 7^ g nr , there is no DFS in the system. In that case, the results 
(not shown) behave similarly to, e.g., Fig. |2j That is, a plot of S ma-X and A± as a function 
of g r , with g r ^ g nr and g p h = 0, shows a clear correlation between the maximum achievable 
value of S and the overlapping resonance contribution. 

V. SUMMARY 

We have introduced a new unifying approach to considering the resistance to decoherence 
of system superposition states. Specifically, using the generic class of spin-boson systems as 
an example we have shown that overlapping resonances of the system bare states, and the 
associated interference between these states, provide insight into the dependence of the deco- 
herence of system superposition states in the important regime of strong system-environment 
coupling. For example, within a given class of Hamiltonians, increasing overlapping reso- 
nances contributions allow for states of higher purity. The approach is completely general, 
although it has been applied here to a spin 1/2 particle interacting with a single thermal 
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oscillator. Extensions to other cases are underway. In the first instance we anticipate that in- 
creasing the environment size will increase the participation and effectiveness of overlapping 
resonances in decoherence control. 
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FIG. 1: Plot of the overlap |-D 7 | 2 between various bare states and the eigenstates I7) as a function 
of the energy £ 7 : (a) | + ,0) , |-,0), (b) | + ,0), |-,2), (c) |+,2), |-,2). The data points for the 
first of each pair are connected by a solid line whereas the data points for the second of each pair 
are connected by solid lines. Data are for the parameters uj = 2ir GHz, loq = 2tt x 100 MHz, 
9r = 9nr = 2ir GHz and g p h = 0. 
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FIG. 2: The dependence of the maximum achievable purity S^ax of the system (solid line) and the 

overlapping resonance Aj (dashed line; right y-axis) on the coupling strength g r /(2ir). Here (a) 

g p h = 2tt x 500 MHz, g nr = 2ir GHz and (b) g p h = 2tt GHz, g nr = 2irx 500 MHz. The remaining 
parameters are u = 2tt GHz, cjo = 2ir x 100 MHz, and the temperature T = 25 mK. The purity is 
calculated at a time t = 0.1 ns. 
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FIG. 3: Variation of purity with time for the initial states (— 0.539998|+)+0.841665|— )) (solid line), 
(1+) + l — ))/v / 2 (dashed line), and |+) (dotted line). Here the interaction strengths g r = 2tt x 400 
MHz, g nr = 2ir GHz, and g p h = 2ir x 500 MHz. The other parameters are as in Fig. [2j 
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FIG. 4: Variation of Sdiff (solid line) and the overlap (dashed line; right y-axis) with g r for 

g nr = 2ir GHz and g p h = 2ir x 500 MHz. The remaining parameters are as in Fig. [2j 




FIG. 5: Variation of S^m (solid line) and the overlap (dashed line; right y-axis) with g r = g ri 

for g p h = 0. The other parameters are as in FigsJ2J 
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